Collapse of the Electron Gas to Two Dimensions in Density Functional Theory 
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Local and semilocal density-functional approximations for the exchange-correlation energy fail 
badly in the zero-thickness limit of a quasi-two-dimensional electron gas, where the density vari- 
ation is rapid almost everywhere. Here we show that a fully nonlocal fifth-rung functional, the 
inhomogeneous Singwi-Tosi-Land-Sjolander (STLS) approach, which employs both occupied and 
unoccupied Kohn-Sham orbitals, recovers the true two-dimensional STLS limit and appears to be 
remarkably accurate for any thickness of the slab (and thus for the dimensional crossover) . We also 
show that this good behavior is only partly due to the use of the full exact exchange energy. 
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One of the main goals of electronic-structure the- 
ory is prediction of the energy of inhomogeneous inter- 
acting many-electron systems. Time-independent and 
time-dependent density functional theories (DFT and 
TDDFT) Q provide calculable predictions respectively 
for ground- and excited-states, and are intimately linked 
together. In these theories, the noninteracting kinetic 
energy is treated as an exact functional of the occupied 
Kohn-Sham (KS) orbitals 0, and only the exchange- 
correlation (xc) energy E xc and/or potential v xc (r) have 
to be approximated. These theories entail a hierarchy of 
approximations for exchange and correlation. The more 
sophisticated approximations satisfy many exact condi- 
tions, but must still be carefully tested for accuracy and 
reliability. In this paper, we test the inhomogeneous 
Singwi-Tosi-Land-Sjolander (STLS) method [3| for the 
quasi two-dimensional (2D) electron gas, a problem that 
becomes increasingly challenging for density functionals 
as the true 2D or zero-thickness limit is approached 

The ladder classification [|| of ground-state density 
functionals for E xc has three complete non-empirical 
rungs: the local-spin-density approximation (LSDA) 0, 
the generalized-gradient approximation (GGA) ||, and 
the meta-GGA |9(. The meta-GGA, which satisfies many 
exact constraints but still uses to some extent the er- 
ror cancellation between exchange and correlation, has 
as ingredients the spin densities n-\ and rij,, their gradi- 
ents Vn-f and Vnj,, and the KS non-interacting kinetic 
energy densities r-f and tj_. These local and semilocal 
density functionals (LSDA, GGA, and meta-GGA) work 
for atoms, molecules, solids, and surfaces 0. They 
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also work for atomic monolayers ll| and other quasi- 
2D systems II 211 . but they fail as the true 2D limit is 
approached [13, EH) because of the high inhomogene- 
ity of the electron density along the confined direction. 
The failure of local and semilocal density functionals to 
describe the dimensional crossover of the exact xc func- 
tional has been avoided by using nonlocal models such as 



the weighted-density approximation 

The next rung of the ladder is the hyper-GGA 
(HGGA) [l6j], a nonlocal- functional approximation which 
uses the Tao-Perdew-Staroverov-Scuseria (TPSS) meta- 
GGA 0] ingredients and the conventional exact exchange 
energies per particle e x -\ and e x i For the xc energy 
per particle at position r, one writes: 
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Here e x (r) is the exact exchange energy per particle in 
the TPSS gauge " 



1 71 ] . and a(r) is a nonlocal functional 
bounded between and 1. The second term on the right- 
hand side of Eq. |T]) is the static correlation and was built 
such that, together with the TPSS (dynamic) correla- 
tion, it is compatible with exact exchange. The mixing 
parameter a(r) goes to 1 when exact exchange is much 
bigger than correlation (for one-electron systems and for 
the high-density limit), when the density is rapidly vary- 
ing, and when an open system has a high fluctuation of 
electron number in spin-polarized regions at the Hartree- 
Fock level. HGGA satisfies more exact constraints than 
any semilocal functional. To balance the full nonlocality 
of exact exchange, universal empirical parameters are in- 
voked in HGGA correlation and fitted to chemical data. 
The set of parameters adopted here is for use with TPSS 
orbitals. We show that this HGGA partially avoids the 
bad behavior of the semi-local xc density functionals for 
the quasi-2D electron gas. 

The central equation of TDDFT linear response (at 
frequency u>), in which all the objects are functionals 
of the ground-state density, is a Dyson-like equation for 
the density-response function x(r, r'; w) (in atomic units 
where e 2 = h = m e = 1) [hil ]: 



X(r, r'; uj) = xo(r, r'; uj) + J dn dr 2 xo(r, ri ; uj) 
— 1 +/sc[n](ri,r 2 ;w)l x(r 2 ,r';a;), (2) 
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where Xo( r , r '; is the density-response function of non- 
interacting KS electrons and is exactly expressible in 
terms of KS orbitals jl9 |. and / xc [n](r,r';u>) is the dy- 
namic xc kernel which must be approximated. When 
f xc [n](r, t';w) is taken to be zero, Eq. © reduces to the 
screening equation of the time-dependent Hartree or ran- 
dom phase approximation (RPA). The xc energy can then 
be calculated using the adiabatic-connection fluctuation- 
dissipation formula [13] , which allows TDDFT to produce 
sophisticated approximations to the ground-state xc en- 
ergy. In particular, the RPA has been evaluated for the 
quasi-2D electron gas [ljl 21 1 and it reaches the RPA 2D 



electron gas limit; however, the RPA xc energy per parti- 
cle e^ c Pj4 for the 2D uniform electron gas underestimates 
the exact e xc by more than 20 mHa/e" [Hj]. In order to 
go beyond the RPA, various approximations have been 
constructed for the xc kernel, but they have mainly been 
taken from the 3D uniform electron gas [23| and do not 
take the dimensional crossover into account. 

The inhomogcneous STLS (ISTLS) 0] is, like the 
RPA it corrects, a fifth-rung density functional that 
employs both occupied and unoccupied Kolm-Sham or- 
bitals. The main idea of ISTLS (and STLS 0) re- 
lies on the truncation of the Bogoliubov-Born-Green- 
Kirkwood-Yvon (BBGKY) hierarchy of the kinetic equa- 
tions [25j, by assuming that the two-particle dynamic 
pair-distribution function (r, p, r', p'; t) can be ex- 
pressed in terms of the single-particle distribution func- 
tion /(r, p; t) as 

/< 2 >(r, p, r', p'; t) = g(r, r') /(r, p; t)J(r', p'; t), (3) 

where g(r, r') is the static (and momentum-independent) 
equilibrium pair-correlation function. 

Using the linearity and time-invariance of the trun- 
cated BBGKY equation, Dobson et al. [3| found the fol- 
lowing Dyson-like "screening" integral equation for the 
density-response function: 



X(r,r';c^) = X o(r,r';w)+ / dr"Q{r, r"; w)x(r", r'; u) 



(4) 



where 



Q(v, r'; w) = - J dr"v (r, r"; w) ■ g(r" , r')V r „ — j- - . 

r (5) 

Here, i>o(r, r'; u) is a vector response function, which sat- 
isfies the equation 



Xo(r,r';w) = V r < • u (r,r';u). 



(6) 



The equilibrium pair-correlation function g(r, r') is ob- 
tained from the fluctuation-dissipation theorem [201 ] : 



fl(r,r') = 1 



7rn(r)n(r') J 



du x(r, r'; iu) 



6{r - r') 
n(r') 

(7) 



Equations (HI)-© are solved selfconsistently, until a con- 
verged solution is obtained. This ISTLS scheme yields 
the exact exchange energy, as does the TDDFT scheme 
described above (see Eq. 0), and only the correlation 
energy is approximated. ISTLS correlation is exact for 
all one-electron densities. 

For the 2D and 3D uniform electron gases, the STLS 
approach made a remarkably accurate prediction of the 
correlation energy over a wide range of densities (for 3D: 
1 < r s < 20 and for 2D: 0.5 < r 2 s D < 16), as confirmed 
by QMC calculations (see Table I of Ref. [26} and ref- 
erences therein) (2^1 ■ The STLS approximation is also 
known to yield reasonable ground-state energies for the 
ID and 2D Hubbard models in the half-filled antifcrro- 
magnetic states j28|. Furthermore, the ISTLS, which is 
a high-level correlated approach that predicts (and does 
not use as input) the correlation energy of the 3D and 2D 
uniform gases, has been shown to yield accurate jcllium 
xc surface energies that are close to their LSDA and RPA 
counterparts, and has been used to demonstrate that a 
local-density approximation for the particle-hole interac- 
tion is adequate to describe the surface energy of simple 
metals [26|. 

A 2D uniform electron gas is described by the 2D 
electron-density parameter r 2 s D — 1/ ^/nn 2D — \/2/ kjP , 
where n 2D is the density of electrons per unit area, and 
k 2 jP represents the magnitude of the corresponding 2D 
Fermi wavevector. The 2D exchange energy per particle 
is 
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A realistic interpolation (which uses QMC data) between 
the high- and low-densit y li mits of the 2D correlation 
energy per particle is [3, H^] 



0.5058 



1.3311 



1 + 1.5026r? D - 1 ) - 



(9) 

For the description of a quasi-2D electron gas, we con- 
sider a quantum well of thickness L in the z-direction. In 
the infinite barrier model (IBM) [30| for a quantum well, 
the KS effective one-electron potential is replaced by zero 
inside infinitely high potential walls at z — and z = L. 
Hence, the normalized KS one-electron wave functions 
and energies at < z < L are 



and 



Ei 



i, fen 



t) 1 



(10) 



(ii) 



where A represents the normalization area in the xy- 
plane, I is the subband index, and r|| and k|| represent 
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the position and the wavevector in the xy plane. In this 
model the electrons cannot leak out of the well, so the 
true 2D electron-gas limit is recovered by simply shrink- 
ing the well. When only the lowest level is occupied 
{E hk 2p < E 2 ,o, i- e., L < y^wr 20 = L max [13]), the 
density of states of this quantum well begins to resem- 
ble the density of states of a 2D electron gas, the motion 
in the z-direction is frozen out, and the system can be 
considered quasi- two-dimensional. 

We would like to contract or expand the electron den- 
sity n(z) without changing the total number of electrons 
per unit area. Hence, we perform a one-dimensional scal- 
ing of the form n\(z) = Xn(Xz) [31], which as A — > oo 
yields the true 2D limit. This scaled electron density co- 
incides with the electron density that one would find from 
Eq. (|10|) by simply replacing the quantum-well thickness 
L by L/X 13]. The corresponding exchange and cor- 
relation energies per particle, £ x [n,\] and s c [n\\ , should 
satisfy the following scaling relations [l3| : 



lim e x [n\] > - 



DC 



> lim £ c [nA] > — oo. 

A — >OQ 



(12) 
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These equations, which start from those of Ref. 
not satisfied by the LSDA, GGA, and meta-GGA. 

Here we investigate the performance of the ISTLS ap- 
proach to describe the IBM quasi-2D electron gas for 
thicknesses L such that L < L max . We choose quasi-2D 
electron gases of fixed 2D electron-density parameters: 
r 2D = 2/\/3 (as in Fig. 1 of Ref. Q1) and r 2D = 4 (as 
in Figs. 2 and 3 of Ref. O). The self-consistent ISTLS 
scheme [Eqs. ((4])-((TJ)] was solved numerically using occu- 
pied and unoccupied orbitals of the form of Eq. (fT0| . 
Our numerical scheme is similar to the one described 
in Ref. [2||; however, in contrast to the jellium surface- 
energy calculation reported in Ref. [2^] , where the ISTLS 
self-consistent scheme converges rapidly, in the case of the 
quasi-2D electron-gas calculation the convergence is slow 
(especially when L << L max ) [32j and special care must 
be taken to perform the frequency integration of Eq. ([7]) . 

Figures [1] and show the results we have obtained, as 
a function of the quantum- well thickness £ m ax/A, for the 
exact exchange energy per particle, s x , and the xc en- 
ergy per particle, e xc , in the following approximations: 
LSDA, PBE-GGA, HGGA, RPA, and ISTLS. We ob- 
serve that, while local and semilocal density-functional 
approximations (LSDA and PBE-GGA) fail badly in the 
zero-thickness (2D) limit, both the RPA and the ISTLS 
nicely recover their 2D counterparts, which in the case 
of the ISTLS (-0.62 Hartree for r 2D = 2/V3 and -0.20 
Hartree for r 2D — 4 [33[ ) ar e very close to the prediction 
of Eqs. ([H])-© or, equivalently, QMC calculations (-0.63 
and -0.21, respectively). The RPA xc energy per parti- 
cle, however, considerably understimates e xc for all slab 
thicknesses 0- Figures CD and [2] show that the ISTLS 
xc energy not only approaches closely the exact 2D limit 
but has also the correct behavior in the limit L — > L max , 
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FIG. 1: Exchange-correlation energy per particle of an IBM 
quasi-2D electron gas of fixed 2D electron density (r 2D = 
2/\/3), as a function of the quantum- well thickness Lmax/A 
(L mal = 4.44). The 2D limit is from Eqs. (©-©. Various 
density-functional approximations have been used (LSDA, 
PBE-GGA, and HGGA), as well as the fifth-rung RPA and 
ISTLS. The exact exchange energy per particle, e x , is also 
plotted, for comparison. While the local and semilocal den- 
sity functionals diverge in the 2D limit, RPA and ISTLS cal- 
culations nicely recover the corresponding xc energy of a 2D 
electron gas. The HGGA recovers the 2D exact exchange 
limit; i.e., as in the case of the LSDA and the PBE-GGA, the 
HGGA correlation energy per particle goes to zero in the 2D 
limit. 



-0.1 
-0.12 
-0.14 
-0.16 
-0.1! 

-0.2 
-0.22 
-0.24 
-0.26 
-0.28 

-0.3 



Exact Exchange 




2D . 

r. =4 



L=L„ 



6 

Jk (bohr) 



8 10 



FIG. 2: As in Fig. [TJ but now for r 2 3 D = 4 (L n 
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which is expected to be well described within the LSDA, 
GGA, and HGGA. The meta-GGA curves, not shown in 
Figs. Q] and [21 are found to be very close to their GGA 
counterparts. We note that the HGGA greatly improves 
over GGA (and meta-GGA) and is more accurate than 
RPA for almost all values of the slab thickness. 

Most of the functionals tested in Figs. Q] and [2] have 
been (exact exchange, LSDA, PBE) or could be (HGGA, 
ISTLS) useful in condensed-matter physics and quantum 
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chemistry. For all of those except LSD A and ISTLS, 
the correlation energy per electron tends to zero as the 
quantum-well thickness £ max /A goes to zero; for ISTLS, 
it tends to a realistic negative value. 

In summary, we have shown that the ISTLS ap- 
proach 1 correctly and remarkably describes the dimen- 
sional crossover (from 3D to 2D) of the xc energy. The 
quasi-2D electron gas is an important and difficult test for 
density- functional approximations. This test is related to 
the one-dimensional scaling, an exact constraint which is 
not satisfied by LSDA, GGA, or meta-GGA. The fourth- 
rung hyper-GGA with full exact exchange (see Eq. (P)) is 
found to improve considerably the behavior of local and 
semilocal density functionals over the whole thickness- 
range of the quasi-2D electron gas. The fifth-rung ISTLS 
scheme, which uses as input all occupied and unoccupied 
KS orbitals and is numerically more expensive [351 ] . is 
found to be remarkably accurate for the description of 
various quasi-2D systems. Due to its accuracy in describ- 
ing both the jellium surface energy [3] and the energy 
per particle of a quasi-2D electron gas, we believe that 
further work in testing and improving the ISTLS scheme 
could be important. So far, the computational cost of 
ISTLS has limited its applications to situations in which 
the density variation is effectively one-dimensional [2h ]. 
Some chemical tests have been carried out, however, at 
the RPA level [36j], and these calculations could well be 
extended within the full ISTLS scheme. 
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